Creating ultracold molecules by collisions with ultracold rare gas atoms in an optical 

trap 
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Abstract 

We study collisions of para-H2 with five rare gas atomic species (He, Ne, Ar, Kr and Xe) over the range from 1 K to 1 /i 
K and evaluate the feasibility of sympathetic cooling Hb with ultracold ground state rare gas atoms co-trapped within a deep 
optical trap. Collision cross-sections over this large temperature range show that all of these species could be used to cool H2 
to ultracold temperatures and that argon and helium are the most promising species for future experiments. 

PACS numbers: 34.50.Cx,37.10.Pq,67.85.-d,34.50.-s 



The ability to trap and cool dilute atomic gases has 
revolutionised physics over the last twenty years. Laser 
cooling has been of primary importance in the explo- 
ration of ultracold interactions between trapped atoms, 
Bose-Einstein condensation, quantum atom optics, and 
for precision metrology. It has also provided a well- 
controlled testing ground for condensed matter physics, 
and more recently quantum information. Considerable 
attention has now turned to the creation and study of 
trapped cold and ultracold molecular gases [lj]. Molecules 
are much more complex than atoms, offering for instance 
the possibility of investigating interactions that are not 
available in atomic species. For example, molecules can 
have a permanent dipole moment that leads to strong 
long-range dipolar interactions, while the anisotropic po- 
larizability of molecules offers the ability to orient them 
in free space and within a trap using electromagnetic 
fields. Cold molecules are also seen as an ideal testing 
ground for physics beyond the standard model of parti- 
cle physics [2| and as a useful tool in the search for parity 
violation at the molecular level. The ability to perform 
chemistry in dilute molecular gases at temperatures far 
below 1 K promises access to chemical processes domi- 
nated by quantum mechanical tunnelling and resonances 
have yet to be observed. 

In general, laser cooling cannot be applied to molecules 
due to the absence of a single cycling transition. Al- 
though ultracold diatomic molecular species can be pro- 
duced by association of laser cooled atomic species on 
Feshbach resonances 0] and by photoassociation 
the range of molecular species that can be produced is 
severely limited by the small subset of atomic species that 
can be laser cooled. New techniques are being developed 
to produce slow cold molecules of much greater variety 
and complexity. Many of these techniques, such as Stark 
deceleration [H, 0], use conservative fields to filter out a 
narrow momentum distribution which can be slowed and 
trapped .7:]. To date all of these approaches been limited 
to temperatures above 1 mK. To produce colder temper- 
atures, without significant loss of molecules, it is nec- 
essary to employ a dissipative cooling scheme once the 
molecules are slowed and trapped. A number of im- 



portant dissipative schemes are potentially available to 
cool molecules and include stochastic, cavity and evap- 
orative and sympathetic cooling. In sympathetic cool- 
ing, the gas to be cooled is co-trapped with a colder gas 
at higher concentration and heat flows from the hotter 
gas to the colder gas, coming to a common temperature 
that is lower than that of the hotter gas. This scheme 
has been demonstrated experimentally using elastic col- 
lisions with cold helium to cool and trap paramagnetic 
molecules in the 100 mK range, with the ultimate tem- 
perature limited by the vapor pressure of helium. More 
recently however, proposals to use this technique to cool 
a range of molecules into the ultracold regime (below 1 
mK) have been explored using ultracold, laser cooled, 
alkali species which can be readily created in magneto- 
optic traps. Theoretical studies to date include collisions 
of OH and NH with laser cooled Rb[l,0]. In these studies 
co-trapping of the ultracold Rb atoms and molecules in 
electrostatic and/or magnetic traps is required. Of crit- 
ical importance to successful sympathetic cooling is the 
dominance of elastic over inelastic processes over a wide 
range of collision energies. Inelastic losses that occur due 
to state changing collisions to untrappable states and/or 
reactions are of particular concern since these processes 
lead to losses from the trap preventing thermalization. 
State changing collisions are particularly important for 
electrostatic and magnetic traps because the molecules 
are not usually trapped in their absolute ground state. 
Finally, the use of an alkali as a collision partner can 
often lead to losses via reactive collisions. For exam- 
ple, Rb-NH interactions can result in a harpoon reaction 
which creates the RbNH molecule [9(. 

We propose a different approach to sympathetic cool- 
ing of molecules that eliminates both of these loss mech- 
anisms and additionally can be used with a range of 
deceleration and trapping methods which have been 
demonstrated to produce stationary, stable and cold 
molecules below 1 K. This method uses inert, rare gas 
atoms as the ultracold collision partner and co-traps the 
molecules and atoms in their absolute ground state in 
a deep optical trap. Ground state rare gas atoms can- 
not be directly laser cooled, but when promoted to their 
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metastable triplet state via electron impact they can 
be laser cooled and trapped in a magneto-optical trap 
(MOT) 0, EE E El ■ Once cooled they can be quenched 
to their ground state via excitation to a higher lying 
state which spontaneously decays to the absolute ground 
state[14|. Cold decelerated molecules can then be spa- 



tially overlapped with the cloud of ground state rare gas 
atoms in a deep quasi-electrostatic optical trap (QUEST) 
or a microwave trap and both species can be held simul- 
taneously. In these traps all states are high field seeking 
so even state changing collisions do not prevent trapping 
and thermalization. However, as relatively warm Stark 
decelerated molecules in the 10-100 mK range must be 
trapped as well, high intensity fields must be used to cre- 
ate a deep trap. 

A buildup cavity offers a route to creating high inten- 
sities and the well depth from such a QUEST is given 
by 
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where the well depth is given by Uq — j^I c , an d a is the 
static polarizability, I c is the one way circulating peak in- 
tensity and k = 2ir/\ is the wavevector along the axis of 
the trap taken to be along x. The 1/e 2 width of the field 
in the radial direction r is given by u and the Rayleigh 
range of the intracavity field is x r = iruj 2 /\. The ra- 
dial and axial oscillation frequencies for particles in the 
trap are given respectively by u r — [4/To/tow 2 ] 1 / 2 and 
lu x = [2C/ fc 2 /to] 1 / 2 . Intensities in the 10 8 W cm -2 range 
can be produced in an optical buildup cavity using laser 
input powers of 17.5 W and a cavity finesse of 6300 [l5j |. 
This intracavity intensity corresponds to a well depth of 
70 mK for C 6 H 6 , 19 mK for NH 3 and 6 mK for H 2 . For 
the ground state rare gases, depths of 27 mK (Xe), 17 mK 
(Kr), 11 mK (Ar), 3 mK (Ne), and 1.4 mK (He) could 
be achieved using this design and it is feasible that much 
deeper optical traps could be created using a near concen- 
tric cavity, but with a reduced trapping volume. The well 
depths achieved for the larger more polarizable molecules 
compare favourably with that achieved to date using elec- 
trostatic traps for molecules and could be combined with 
these traps to increase phase space density. The deeper 
well depths also lead to higher oscillation frequencies 
which will reduce the time to thermalize species withi 
the traps. As rare gas atoms in their ground state have a 
relatively high ionization potential (lowest is 12.1 eV for 
Xe) they are robust against multiphoton/tunnel ioniza- 
tion within the high intensity fields of the trap. Similarly, 
most stable molecules of interest also have high ioniza- 
tion potentials and no ionization is seen using intensities 
even in the 10 12 Wcm~ 2 range The choice of which 
rare gas species to use for sympathetic cooling must de- 
pend primarily on favorable elastic to inelastic collision 
cross sections between the species, but also on the po- 
larizability of both the atom and molecular species. Of 



the rare gas atoms considered, xenon is the most polariz- 
able and the deepest traps can be created for this species. 
However, the lighter rare gases such as helium and neon 
may be more favorable for cooling lighter molecules due 
to the favourable mass ratio and weaker Van der Waals 
interactions. 

The aim is develop a method for trapping arbitrary, 
stable (ie closed shell) molecules. To explore the feasi- 
bility of this scheme we have chosen to study collisions 
of H2 with ultracold rare gases. Study of these collisions 
is important for their fundamental nature, but also for 
their importance in testing ab initio theoretical models 
and their relevance for astrophysics. H2 is also seen as a 
ideal candidate to study chemistry at ultracold tempera- 
tures dominated by resonance and tunneling phenomena 
dEEl- Cold H 2 in the 10 mK to 100 mK range can in 
principle be decelerated using optical Stark deceleration 
|16l ll9[ from an intense molecular beam which can be 
created in pure para form with initial temperatures in 
the 100 mK range. 

Cold and ultracold collisions and calculations of col- 
lision cross sections for H 2 with He and Ar have been 
already considered in some detail 

USUI- In 

one sense, 

we extend these to include all stable rare gas atoms that 
can be laser cooled and trapped. On the other hand, 
much of the previous work focused on internally highly 
excited H 2 molecules, and was directed to study rovibra- 
tional energy transfers. For this purpose, a system of 
coupled equation needs to be solved containing at least 
as many equations as the number of open channels 21 1. 
The experimental setup proposed aims at producing H 2 
that is decelerated for example by optical Stark deceler- 
ation in its vibrational-rotational ground state. Under 
this very specific experimental condition, the scattering 
of Rg-H 2 can be effectively be modelled as a two-body 
process [22j . As the scattering energy is very small com- 
pared to the vibrational-rotational excitation energies of 
H 2 , the only allowed channel is the elastic one, and no 
inelastic processes need to be considered. Under those 
assumptions, the Rg-H 2 system can be effectively mod- 
elled as a two-body system 24|, with an effective Rg-H 2 
interaction 



V eff (R) = J 2 o (r)P 2 ( M )U Rg _ H2 (r,i?,M)^, (2) 

where r is the H-H distance, R is the Rg-H 2 dissociation 
coordinate, fx = f ■ R and Pq is a Legendre polynomial. 
</>oo represents the ground vibrational-rotational function 
of H 2 , and VR g _H 2 ( r : R> A 4 ) is the R-gEb potential energy 
surface. From the asymptotic behaviour of the Rg-H 2 
wavefunction at large R one can extract the phase-shift 
8(E), thus obtaining the elastic cross-section a(E). The 
zero energy limit of cr(E) is well known, and for s— wave 
scattering is a(E = 0) = 47ra 2 , where a s is the Rg-H 2 
scattering length. Therefore the determination of the 
scattering length provides a functional tool to compare 
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the ultracold scattering properties of the different com- 
plexes. 

The scattering length can be dramatically sensitive to 
small variations in the potential energy surface (PES). 
This is particularly true in systems which are very weakly 
bound. Therefore, the accuracy of the underlying molec- 
ular PES is essential for a sensible determination of the 
Rg-H2 a s . As a consequence of a long standing inter- 
est in Rg-H2 complexes, the Rg-H2 interaction has been 
studied in great detail, in particular for ArH2 and HcH2 . 
However, the Rg-H2 Van der Waals interaction is weak 
due the closed shell structure of both H 2 and the rare gas 
atoms. Despite many years of work, there is still a small 
but critical uncertainty associated with the Rg-H 2 PES's. 
For example, in HeH2, which is the weakest bound of all 
five complexes, the existence of a bound state has been 
long debated, and only recently the calculated potentials 
have become precise enough to confirm its existence [23| . 
To estimate the scattering lengths of the various com- 
plexes, we have gathered several available PESs from the 
literature. Table [J reports the range of values obtained 
for the Rg-H2 scattering length and effective range. Of 
all complexes and isotopologues, 3 HeH 2 represents a spe- 
cial case as the percentage variation in its mass makes its 
collisional properties sufficiently different from 4 HeH 2 to 
make it worth mentioning explicitly. The values of ta- 
ble [J have been calculated using a single channel trial 
wavefunction, and corrections were also calculated and 
proved insignificant 
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Balakrishnan et c 
performed a detailed calculation of the He-Eb collisional 
process using both PES's of Ref. and Ref. [H; 
our results agree well with theirs. Table |l] shows that, 
although there are big uncertainties in the collisional pa- 
rameters associated with the uncertainties in the PES, 
the physics of the scattering mechanisms is rather clear. 
The greatest cross-section is of He-H 2 , due to the strong 

' US]. The con- 



halo characteristics of this system [22 
cept of quantum halo systems is borrowed from nuclear 
physics [29( and defines quantum systems which extend 
well into the classically forbidden region. It is possible 
to show that for those systems the scattering length de- 
pends on the square root of the inverse of the binding 
energy E, a s w l/\/E, therefore the smaller the binding 
energy the larger the scattering length. The scattering 
lengths of the five complexes can be interpreted in terms 
of this relationship. For example, the two HeH 2 isotopo- 
logues are by far the weakest bound of the series thus 
their scattering length is the greatest. The NeH 2 system 
is more strongly bound, therefore its scattering length is 
smaller. The stabilization of a second vibrational band 
in ArH2 makes this complex more reactive than its pre- 
decessor. Again, as the potential binding increases, the 
KrH2 a s is smaller than the one of ArH2. The Xe-H2 
scattering length is the smallest of the whole group. 
The XeH2 PES is the one with the greatest uncertainty. 
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FIG. 1: Collisional cross-sections for Rg-H 2 complexes over 
the energy range from 10 K to lfi K 
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re// [A] 


a(E = 0) [A 2 ] 


PES Refs. 


3 He-H 2 


67.6-90.6 


6.08-6.27 


57500-103000 


[27, 28, 30] 


4 He-H 2 


22.7-24.7 


5.00-5.16 


6500-7800 


[27, 28, 30] 


Ne-H 2 


3.30-3.85 


10.7-14.9 


140-190 


[30, 31, 32] 


Ar-H 2 


8.71-10.1 


5.21-6.12 


950-1300 


[30, 33, 34, 35] 


Kr-H 2 


5.51-6.96 


9.65-14.8 


380-610 


[33, 36, 3JJ 


Xe-H 2 


1.82 


289 


42 


f33] 



Tests on the sensitivity of the results of Table U to the 



TABLE I: Ultracold collisional properties for Rg-H 2 com- 
plexes, namely the scattering length a a and the effective range 
r e //. The zero-energy elastic cross section a is also tabulated. 



potential showed that Ne-H2 and KX-H2 are essentially 
unaffected by increases or decreases of the interaction 
potential by over 25 %. However the a s for Ar-H2 and 
Xe-H2 can become comparable to the one of 4 He-H2 with 
relatively small adjustment of their PES, —25 % and + 
25 % respectively. As the Ar-H2 PES is weakened, the 
already weakly bound excited vibrational band becomes 
closer to the dissociation threshold, whereas for the Xe- 
H2 strengthening the potential brings a third vibrational 
state close to being bound. 

The large collisional cross-sections for 4 He-H2 and Ar- 
H2 are comparable to that utilized in sympathetic cooling 
experiments between cold alkalis [38] and these values are 
relatively constant over a large temperature range. We 
estimate a collision rate and thermalisation time for Ar- 
H2 based on a cross-section of 1000 A 2 and radial and 
axial trap frequencies of 44 2ir kHz and 37.4 2ir MHz 
for Ar and 16.5 2tt kHz and 13.8 2tt MHz for H 2 . We 
also assume that H2 has evaporated to approximately 
1/10 of the well depth before sympathetic cooling oc- 
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curs. We assume that each species has a density in the 
trap of the form n = no exp [— U/kT] with peak densi- 
ties of 5 x 10 9 cm -3 and 10 11 cm~ 3 corresponding to 11 
molecules and 1300 atoms per fringe of the trap at initial 
temperatures of 600 /iK and 10/^K respectively. The Rg- 
H2 collision rate is given by 7 = av re l J nR g (x)rtH 2 (x)dx, 
v re i is the mean relative velocity between collision part- 
ners and riRg (x) and nn 2 ( x ) are the densities of the rare 
gas and Hi- It is well known that approximately 3 colli- 
sions are required for thermalization of collision partners 
of equal mass in a gas but for unequal masses this is ap- 
proximated by - where ri — 4 m Rg m H 2 rp^ initial 

time for thermalization for two species within a trap is 
given by [IH 

t 37r 2 fc b r H2 , 3 « 

For the conditions given for Ar-H 2 , a thermalization time 
scale of 300 ms is determined. This relatively fast ther- 
malization time can be accounted for by the large axial 
trap frequency and the large relative velocity (2.5 ms -1 ) 
between the hot molecular gas and the much colder 
atoms, despite the poor spatial overlap between the two 
species. For Xe-H2 collisions, using the same initial den- 
sities and a collision cross-section of 42 A 2 , a thermal- 
ization time of approximately 30 seconds is estimated. 
Although the Xe-H2 collision rate is much lower than for 
Ar-H2 trapping of atomic species over this time scale has 
been demonstrated indicating that even Xe, which forms 
the deepest trap, could be used for sympathic cooling. 

We present a new route towards the creation of ultra- 
cold molecules by sympathetic cooling with laser cooled 
metastable rare gas atoms that are quenched and subse- 
quently co-trapped in a deep optical lattice trap. The 
scheme allows the simultaneous trapping of relatively 
'warm' molecules and ultracold rare gas atoms and uti- 
lizes the high trap frequencies in deep standing wave 
traps for increasing the collision rate and reducing ther- 
malization times. In order to evaluate this scheme for 
cooling molecular hydrogen we have calculated the colli- 
sional properties for all five rare atoms. Using the cross- 
sections determined in this study we show that sympa- 
thetic cooling of H2 molecules by all five atoms is tech- 
nically feasible over a wide temperature range from the 
cold to ultracold regime and that this method may be 
promising for a large range of molecular species. 
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